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Abstract. Let A2{t) be the Airy2 process. We show that the random variable 

sup{A(i) -t^} + mm{0,Q}^ 



I has the same distribution as the one-point marginal of the Airy2-).i process at time a. 

■ These marginals form a family of distributions crossing over from the GUE Tracy- Widom 
^ I distribution -FGUE(a;) for the Gaussian Unitary Ensemble of random matrices, to a rescaled 

, version of the GOE Tracy- Widom distribution -FGOE(4^'^^a;) for the Gaussian Orthogonal 

■ Ensemble. Furthermore, we show that for every a the distribution has the same right tail 

■ decay e~'s^^ . 



1. Introduction 



Oh 

■ The Airy processes are a collection of stochastic processes which are expected to govern 

. the spatial fluctuations of random growth models in the one dimensional KPZ universality 

class for wide classes of initial data. They are defined through their finite dimensional dis- 
tributions, which are given by Fredholm determinants. The three basic processes are Airy2 
[PS02], corresponding to curved, or droplet initial data; Airyi [Sas05; BFPS07; BFP07], 
^> I corresponding to flat initial data; and Airystat [BFPIO], corresponding to equilibrium initial 

^ ■ data. 

The KPZ class is identified at the roughest level by the unusual t^^'^ scale of fluctua- 
CN I tions. It is expected to contain a large class of random growth processes, including the 

Kardar-Parisi-Zhang equation itself, as well as randomly stirred one dimensional fluids, 
polymer chains directed in one dimension and fluctuating transversally in the other due to 
a random potential (with applications to domain interfaces in disordered crystals), driven 
lattice gas models, reaction-diffusion models in two-dimensional random media (including 
biological models such as bacterial colonies), randomly forced Hamilton- Jacobi equations, 
^ ■ etc. A combination of non-rigorous methods (renormalization, mode-coupling, replicas) 

and mathematical breakthroughs on a few special models has led to very precise predic- 
tions of universal scaling exponents and exact statistical distributions describing the long 
time properties. These predictions have been repeatedly confirmed through Monte-Carlo 
simulation as well as experiments; in particular, recent spectacular experiments on turbu- 
lent liquid crystals by Takeuchi and Sano [TSIO; TS12] have been able to even confirm 
some of the predicted fluctuation statistics. 

The conjectural picture that has developed is that the universality class is divided into 
subuniversality classes which depend on the initial data class, but not on other details of the 
particular models. Because of their self-similarity properties, the three basic initial data 
are, at the level of continuum partition functions (taking logarithms gives free energies 
or height functions): Dirac ^Oi corresponding to curved, or droplet type initial data; 0, 
corresponding to growth off a fiat substrate; and e^^^^ where B(x) is a two sided Brownian 
motion, corresponding to growth in equilibrium. Of course, in discrete models of various 
types one is dealing with discrete approximations of such initial data. There are also 
three additional non-homogeneous subuniversality classes corresponding to starting with 
one of the basic three on one side of the origin, and another on the other side. The spatial 

1 



SUPREMUM OF THE AIRY2 PROCESS MINUS A PARABOLA ON A HALF LINE 



2 



fluctuations in these six basic classes of initial data are supposed to be given asymptotically 
by the six known Airy processes: the three basic Airy processes, Airy2, Airyi and Airystat) 
and the crossover Airy processes Airy2-i>i [BFS08], Airy2_5.BM [IS04; CFPIO] and Airyi_5.BM 
[BFS09]. 

However, since all initial data are superpositions of Dirac masses, there is a sense is 
which the Airy2 process is the most basic. Although the various microscopic models are 
not linear in the initial data, this is the case for the stochastic heat equation, whose 
logarithm is the solution of the KPZ equation. And for other models, the linearity should 
hold asymptotically. In the limit, the logarithm of the superpositions of exponentials of 
Airy2 processes becomes a variational problem. 

The conclusion is a conjecture that the one-point marginals of the other Airy processes 
should be obtained through certain variational problems involving the Airy2 process. The 
first example of this was the celebrated result of Johansson [Joh03] (see also [CQR12]) that 
the supremum of the Airy2 process minus a parabola has the same distribution as a rescaled 
version of the one-dimensional marginal of the Airyi process, i.e., the GOE Tracy- Widom 
distribution: 

(1.1) P( sup {A2{x) -x^} <m] = P(^i(0) < 2-^/'^m) = FcoEi^^^^rn). 

The general conjecture is based on heuristics which we describe next in the context of the 
stochastic heat equation. 

1.1. Heuristics. We will explain the heuristics first for the case of the Airy2-j.i process. 
Let z(t, x) denote the solution of the one-dimensional stochastic heat equation 

dtz = \d1z — 

where .^(t, x) is space-time white noise. The solution at position x and time t starting from 
a Dirac mass at y at time can be written as 

(1.2) .(0,2/;t,x) = _^e-T-^+2-W/3^,(2-i/3,-./3(._,))^ 

where At is conjectured to converge to the Airy2 process, At{x) — )• A2{x) (see Conjecture 
1.5 in [ACQll] for a precise statement, and [PSll] for a non-rigorous derivation). Starting 
from the step initial data z{Q, x) = lx>o the prediction is 

(1.3) -\ogz{t,x) ^ ^x2l^<o + ^i + log(^/2^)-2-W/M2^i(2^^/V2/3x). 
On the other hand, by linearity we have for each fixed x, in distribution, 

(1.4) .(t, X) = dy .(0, t, x) = dy _i^,-^-^.2-V3,V3^,(2- V3,-./3(._,)) _ 

Note however that as written the equality can only hold in distribution for each t and 
X. If one wants a stronger statement, for fixed t but multiple x, one has to replace 
^t(2"^/^t"2/3(x - y)) in (1.2) by a two parameter process At{2-'^/^t~'^/^x,2-'^/^t~'^/^y), 
keeping track of the statistical dependence on the initial y. For fixed y, it is distribution- 
ally, as a process in x, equal to At{2~^^^t~'^^^(x — y))- And, by symmetry, the same is true 
for fixed x, as a process in y. However, they are not equal in distribution in the sense of 
two parameter processes in both x and y. The limit of At{x,y) is unknown at this time, 
so one is stuck at the level of one-dimensional marginals. 

Calling X = 2~^/^t~^/^x and y = 2~^^^t~'^^^y we can rewrite the exponent in (1.4) as 

2-W/3[^,(x-y)-(5-y)2] 

so that for large t the fiuctuation field 2^/'^t^^/^ [log z{t,x) + -^t + log(\/2vrt)] is well 
approximated by 

sup {A2{x -y) - {x - yf). 

y>o 
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Comparing with (1.3) we deduce that the processes supy>o (^2(2; — y) — {x — y)^) and 
•^2~>i{x) — x^lx<o should have the same one-dimensional distribution or, equivalently, that 

(1.5) A2^i{x) - x^lx<o = sup{^2(y) - y^} 

y<x 

for each fixed a; € M. 

The same argument works for the other two crossover cases. If we let z(0, x) = e^'^^'>\x>o-: 
where B{x) is a standard Brownian motion, then (1.3) and (1.4) are replaced respectively 
by 

- log z{t, x) ^ ^^xH.^o + jgt + log(V2^) - 2-i/3ii/3^2_^j^j^(2-i/3i-2/3^^ 

and 

zit,x) = rdyz{0,y;t,x) = f i ,--^-^.b(.).2-V3,v3^,(,-V3,-./3(._,,^ 



J u 

and now the same scaling argument allows to conjecture that 

^2^bm(2;) - x^la;<o = sup (^2(y) + B{x - y) - y^) 

y<x 

for each fixed rc G M, where now B{y) is a Brownian motion with difi'usion coefficient 2. 
An analogous argument with z{0,x) = lx<o + e^^^^lx->o translates into conjecturing that 

^i^BM(a^) = sup {A2{y) + B{x - y)ly<x - y"^) 



for each fixed rr S M. 

In this article we prove the conjecture (1.5) for ^2-5.11 which connects the three Airy 
processes with non-random initial data. To state the result precisely, we now recall the 
exact definitions of the Airy2, Airyi, and Airy2_j.i processes, together with some additional 
background. 

1.2. Statement of main results. The Airy2 process A2, introduced by Prahofer and 
Spohn [PS02], is a stationary process on the real line whose one dimensional marginals 
are given by the Tracy- Widom largest eigenvalue distribution for the Gaussian Unitary 
Ensemble (GUE) from random matrix theory [TW94]. It is expected to govern the asymp- 
totic spatial fluctuations in a wide variety of random growth models on a one dimensional 
substrate with curved initial conditions, and the point-to-point free energies of directed ran- 
dom polymers in 1-1-1 dimensions (the KPZ universality class). It also arises as the scaling 
limit of the top eigenvalue in Dyson's Brownian motion [Dys62] for GUE (see [AGZIO] for 
more details). It is deflned through its finite-dimensional distributions, which are given by 
a determinantal formula: given ^1, . . . , G and ti < • • • < in M, 

nMtl) < 6, • • • ,A2{tm) < U) = det(/ - f'/'i^2"*f'/')i2(|4^,...,4„}^K)' 

where P denotes probability, we have counting measure on {ti, . . . , t^} and Lebesgue mea- 
sure on M, f is defined on {ti, . . . , t^} x 1^ by 

(1-6) f(tj,x) = la.g(^^.^oo), 

and the extended Airy kernel^ [FNII99; Mac94; PS02] is defined by 

f dX e-^(*-*') Ai(e + A) Ai(e' + A) , \it>t' 
\- /\ dAe-^(*-*') Ai(^ + A) Ai(e' + A), if t < t', 

where Ai(-) is the Airy function. 

The Airyi process, introduced by Sasamoto [Sas05], is another stationary process, whose 
one-point distribution is now given by the Tracy- Widom largest eigenvalue distribution 
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for the Gaussian Orthogonal Ensemble (GOE) from random matrix theory [TW96]. It is 
defined through its finite-dimensional distributions, 

nMh) < ei, . . . ,^i(tn) < Cn) = det(/ - f^/'i^rf^/')^.„,^,...,„},M), 

with f as in (1.6) and 



^2 



'>t 



\/47r(t' - t) "V 4(t'-t) 

+ Ai(e + e + it' - tf) exp((i' -m + i') + ^(t' - tf^ . 

It is expected to govern the asymptotic spatial fluctuations in random growth models with 
flat initial conditions, and the point-to-line free energies of directed random polymers. 

The Airy2^.i process .42^.1, introduced by Borodin, Ferrari, and Sasamoto [BFS08], is 
given by 

nM^lih) < el,...,^2^l(^„^) <em) =det(/-f^/2^oof'/')^2(|,^,,„,,,,^|^K)' 

with f as in (1.6) and 

1 f (y - 

Koois,x;t,y) = , , exp -— lt>, 

V47r(t - s) V 4(t -s)J 

+ — rTTT / dw / dz ■ 



(27ri)2 J^_ gz^/s+sz^-xz -w){z + w)' 

where x = x — s^ls<0; V = U — t'^^t<o a^id the paths 7+, 7- satisfy —7+ C 7„ with 
7+ : e''^+oo e-'"t'+oo, 7_ : e'^'^-oo e''^-oo for some (p+ G (7r/3,7r/2), £ (7r/2,7r - 
(/)+). The Airy2_5.i process crosses over between the Airy2 and the Airyi processes in the 
sense that A2->i{t + T) converges to 2^/^^i(2^^/^r) as t ^ 00 and ~4,2(t) when t — —00. It 
is expected to govern the asymptotic spatial fluctuations in random growth models when 
the initial conditions are half flat. In particular, it is shown in [BFS08] that it governs the 
asymptotic fluctuations for the totally asymmetric (to the left) simple exclusion process 
starting with particles only at the even positive integers. 

Define the crossover distributions C^^, for a G M, as follows: 
(1.7) Gl-^\m) = det(/ - PmKaPm) 

where Pm denotes the projection onto the interval [m, 00), Ka = + and the kernels 
and are given by 

/•oo 

Kl{x,y)= / dAe2"-^Ai(x- A-Fmax{0,a}2)Ai(y-h A-hmax{0,a}2) 
JO 



and 



Kl{x,y) = / (iAAi(x-F A-Fmax{0, a}^)Ai(y + A-hmax{0,a}^). 
Jo 



Here, and in everything that follows, the determinant means the Fredholm determinant 
in the Hilbert space L^(M). As noted in Appendix A of [BFS08], the kernel K^o can be 
expressed in terms of Airy functions^: 

(1.8) K^is, t; X, y) = Lo{s, x; t, y) + e2*V3-2s3/3+ty-sx[^^ ^ ^. 



This corresponds to a minor correction of the formula appearing in [BFS08], where the exponential 
prefactor appears in front of Lq instead of Li +1/2- 
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where 

Lois,x;t,y) = -e(-*)^(x,y) = - ^ -i^-yfm-s)^ 

V47r(t - s) 

/•oo 

Li{s,x;t,y) = / dA e^('+*) Ai(x - A) Ai(y + A), 
Jo 

/•oo 

L2{s,x,t,y)= / (iAe^(*-^)Ai(x + A)Ai(y + A) 
Jo 

with X = X - s^ls<o, y = y - t^lt<o, X = X + s^ls>o and y = y + t^lt>o- Using this for 
s = i = a it is straightforward to check that Koo{t, -'yt, ■) is just a similarity transform of 
the kernel Ka, and therefore 

Gl-^\m) =F{A2^iici) < m). 

Remark 1.1. The operator PmKaPm appearing inside the determinant defining G'^^ in 
(1.7) is trace class (this follows from (2.12) together with a similar bound for K'^). This 
should be compared with the fact that the extended kernel given in [BFS08] for the higher 
dimensional joint distributions of ^2->i is not trace class (but, as shown in Appendix B of 
[BFS08], there is a conjugate kernel which is). 

The following result confirms the conjecture (1.5): 
Theorem 1. Fix a G M. For every m G M, 

f(^sup {A2{t) -f)<m- min{0, q}^^ = G'^^^{m). 

We remark that the equality is easy to obtain in the limits q — t- 00 and a — ^ —00. 
Since the Airy2 process is stationary one expects that, as a — )• —00, sup^<(;j (-^2(0 — t^) is 
attained at t ~ a, and thus supj<Q, (^2(i)~^^) +niin{0, a}^ ~ ^2(0), which has distribution 
-^GUE("i). For the right hand side, observe that for a < the kernel K"^ equals the Airy 
kernel Kx\, which is given explicitly by 



/•oo 

:i.9) KAi{x,y)= / (iAAi(x + A)Ai(y + A). 

Jo 



In addition, one can check that, due to the exponential factor e^"'*', P^K^Pm goes to in 
trace norm as a — —00 (see (2.16)). Consequently, 

(1.10) lim G^-^i(m) = det(l-P^KAiPm) 

a— >— cxD 

which is also Fgue("i). 

On the other hand, as a — )■ 00, the left hand side becomes P(supjgi[j {A2{t) — t^) < m) = 
FGOE(4i/3m). For the right hand side, it is not hard to check that PmK^Pm goes to in 
trace norm as a — )■ 00. In addition, using (1.8) and (A. 6) of [BFS08] with ri = T2 = a we 
get for a > 

K^ix, y) = 2-1/3 Ai(2-i/3(x + y)) - K^x, y) 

with K^{x,y) = dXe^'^^ Ai(x+a^ + A) Ai(y+a^ — A), and one can check that PmK^Pra 
goes to in trace norm as a — )■ 00 (see the comment following (2.16)). The first term on 
the right hand side above corresponds to the kernel B(x, y) = 2~^^^ Ai{2~^^^{x + y)), and 
we deduce after changing variables x 1— )• 2~'^^^x, y 1— )• 2~'^^^y in the resulting determinant 
that 

(1.11) hm Gl-^\m) = det{l - P^^/3„,BP^^f^J , 

where B{x, y) = ^ Ai(i(x + y)), which is also equal to FQOEi^^^^'m) [FS05]. 

The fact that G"^^ crosses over between the GUE and GOE distributions is of course 
a particular case of the crossover property of the Airy2_s.i process. Note that the scaling 
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by 4^/^ in the GOE end of this interpolation implies that both ends satisfy the same 
asymptotics log (l — G'^^{m)) ~ — |m^/^ as m — 00^. In fact, the same upper bound 
holds for ah a G M: 

Proposition 1.2. For every a £ M, there is a c> such that, 

< cme"i™'^'+"™ asm ^00. 

The proof of Theorem 1 is based on a continuum statistics formula for the Airy2 process, 
developed in [CQR12], which is well adapted to such variational problems. Fix a function 
g G H^{[£,r]) and introduce an operator 0^^^] which acts on L^(M) as follows: ©^^^j/C') = 
u{r, •), where u{r, •) is the solution at time r of the boundary value problem 

dtu + Hu = for X < g{t), t E (£, r) 

(1.12) u{a,x) = f{x)l^^g(^a) 

u{t, x) = for X > g{t) 

for the Airy Hamiltonian, 

H = -dl + X. 

The formula reads 

(1.13) nMt) < git) for t G [i,r]) = det(/ - i^Ai + ef^^^f^'-'^'' Km 

Choosing g{t) = m — min{0,a}^ + and r = a, (1.13) gives an explicit formula for the 
probability in Theorem 1 in the limit ^ — )• — 00. The proof will consist on computing this 
limit and showing that it coincides with G^^^(m). 

Note that this strategy is considerably more difficult to implement for the Airystati 
Airy2^BM and Airyi^BM processes, because it involves computing an expectation of the 
Predholm determinant in (1.13) with respect to Brownian motion paths. For example, in 
the stationary case, for which the one-point distribution is given by ¥{Asta.t{t) < '^) = 
i<o(m) with Fq the Baik-Rains distribution [BROO], we would need to check the formula 

(1.14) Fo(m) = hm Efdetfl - i^Ai + Of^'l^^^^'^^'e^^^ Km 



L— >oo 



where 0[_y 2,] explicit kernel, which can be derived from Theorem 3 of [CQR12] , 

and is given by 

(1.15) 

3/)-^ /8L 

eP^^tr^'"(^'?/) = e-^(-+^)+^^V3 ^ _ p-^,^,^^,^, (Bis) < B{s) + m on [-L,L] 

where B a Brownian bridge from x — L? at time —L toy — L? at time L, B is an independent 
two sided Brownian motion with -B(O) = (both with diffusion coefficient 2), and the 
expectation E in (1.14) is with respect to B while the probability P in (1.15) is with 
respect to B. 

1.3. Connection with last passage percolation. It is worth remarking that the general 
picture we have described holds essentially exactly at the discrete level in the case of last 
passage percolation. Here one considers a family {w{i, j)}ij(zi+ of independent identically 
distributed random variables and lets n„ be the collection of up-right paths of length n, 
that is, paths tt = (vro, . . . , 7r„) such that tTj — -ki-i G {(1,0), (0, 1)}. The point-to-point last 
passage time is defined, for m,n £ Z^, by 

m+n 

LP°'"Vm,n)= max wMi)), 

7ren,„+„:(0,0)->{m,n) ^ 



^This is due to the known asymptotics log(l — -Fgoe("i)) ~ — |m''''^ and log(l — FGUEim)) ~ —^m'^^^, 
which follow from the formulas for these distributions in terms of the Painleve II function [TW94; TW96]. 
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where the notation in the subscript in the maximum means ah up-right paths connecting 
the origin to {m,n). Similarly, the point-to-line last passage time is defined by 

L^'^'^n) = max LP°'"*(^ -k,n + k). 

fc=— n,...,n 

Next one defines the process 1 1— )■ i?n°™*(t) by linearly interpolating the values given by 
scaling LP°'"*(m, n) through the relation 

^point^^ + y,n - y) = cin + C2n^/3i/P™"t(c3n^2/3y)^ 

where the constants Cj depend only on the distribution of the w{i,j). In the special case 
where the w{i,j) have a geometric distribution, Johansson [Joh03] showed that 

(1.16) HP^'^^'it) ^ A2{t) - t^ 

in distribution, in the topology of uniform convergence on compact sets, where A2 is the 
Airy2 process. One can also define a rescaled version H^^'^ of L^™'^(n), obtaining the relation 

i^«°« = sup{i^P™°*(t)} 

(here we are setting ffn°™*(i) = for \t\ > csn^/^). It is known [BROl] that H^l'^^ converges 
in distribution to a GOE Tracy- Widom random variable, and hence (1.16) allows to take 
n — )• 00 in the last equality to recover (1.1) (this was Johansson's original proof of (1.1), 
the proof in [CQR12] is based on (1.12)). 

In principle, this idea can be extended to the obtain variational formulas for the other 
Airy processes. For example, one could attempt to replace point-to-line last passage times 
by point-to-half-line last passage times to recover (1.5). Unfortunately, the connection 
between Airy2_j.i (and the other Airy processes) and last passage percolation is made 
through translating the corresponding results for the totally asymmetric exclusion process, 
and by doing this the boundary conditions end up away from the line {(n — k,n -\- k), k = 
—n, . . . , n}, so (1.16) is not directly applicable. In work in progress Corwin, Liu, and Wang 
[CLW] obtain an improved version of the slow decorrelation result proved in [CFP12], which 
would lead to a general version of formulas for last passage times in last passage percolation 
in terms of variational problems for the Airy2 process. In particular, such a result would 
give a proof of the conjectures made in Section 1.1 and, as a consequence, would show that 
(1.14), and similar formulas for the other Airy processes, hold. 

Acknowledgements. JQ and DR were supported by the Natural Science and Engineering 
Research Council of Canada, and DR was supported by a Fields-Ontario Postdoctoral 
Fellowship and by Fondecyt Grant 1120309. Part of this work was done during the Fields 
Institute program "Dynamics and Transport in Disordered Systems" and the authors would 
like to thank the Fields Institute for its hospitality. 

2. Derivation of the formula 

As in [CQR12; MQR12] we will first give an expression for the distribution of the supre- 
mum over a finite interval and then take a limit. Thus we choose an L > —a, which will 
later be taken to infinity, and work on the interval [— L,a]. For notational simplicity we 
will write 

m = m — min{0, a}^. 

We recall that the shifted Airy functions (t)x{x) = Ai(x — A) are the generalized eigenfunc- 
tions of the Airy Hamiltonian, as H(f)x = X4>\, and the Airy kernel K^i is the projection 
of H onto its negative generalized eigenspace (see Remark 1.1 of [CQR12]). Therefore 
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g(a+L)H^^. j^^g integral kernel 

(2.1) e^^+^^^KAiix^y) = / dXe^^^^^^ Ai{x - X) A[{y - \). 

J —00 

We also deduce that e^''+^^^ Km = KAie^'^^^^^ Km, so we may use the cyclic property of 
determinants to rewrite (1-13) as 

(2.2) F{A2{t) < g{t) for t G q]) = det(/ - Km + e("+^)^AAief_^^^jKAi) • 

We will apply the above for g{t) = + m, and for this choice of g we will write @a,L for 
®f-L a]' case the resulting kernel can be computed explicitly, and a minor variation 

of (1.4) in [CQR12] gives 

where Pm = I — Pm denotes the projection onto the interval (— oo,m] and Ra,L is given by 

(2.3) RaL(x,y) = , e ^ ^ " i(o,+L) 
^ ^ ^ V47r(a + L) 

Following [CQR12] we decompose Qa,L as 

R , — p-{'^+L)H p „ _ /? .p O r 

where ^la,L = Pm+L^{^~^"^^^^ ~ Poi,L)Pm+a^- Using this in (2.2) we can write 

(2.4) F{A2{t) < for t G [-L, a]) 

= det(/ - KMPm+a^KM - 6^''+''^'' KMRa,LPm+a^KM - e("+^)^irAi^^a,LAAi) . 
We will show below that 

(2.5) na,L := e("+^)^i^AiJ^a,Li^Ai > 

L— >oo 

in trace norm. Then since the mapping A 1— )• det(/ + A) is continuous in the space of trace 
class operators (see (2.11) below), all that is left to do in order to take L — )• 00 in (2.4) is 
to compute the limit of the operator e^°'~^^^^ KAiRa,L 

We will first proceed formally to identify the limit, and then verify it in Lemma 2.1. 
Since Kai is a projection and H leaves Ka\ invariant, we will pretend that and 
Ka\ commute, so we have to compute the limit of 6^°^^^^^ Ra^i- Define the reflection 
operator Qm by 

Qnif{x) = /(2m - x). 
Then the operator Ra,L defined in (2.3) can be rewritten as 

(2.6) R^,L = e("'+^')/=^e-^«e(-+^)^^>^+(^.+,.)/2e-«. 

Here e^^ (^ stands for a generic variable) denotes the multiplication operator (e^'^/)(x) = 
e^^ f{x) and A is the Laplacian, so e^^ is the heat kernel as in the introduction. The 
advantage of this form for Ra,L is that it will allow us to use the Baker-Campbell-Hausdorff 
formula to compute formally e^°'~^^^^ R^^l- 

We will use the following identities, where [•, •] denotes commutator: 

[H, A] = A] = -2V, [H, V] = [e, V] = -/, [H, i\ = -2V. 
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If A and B are two operators such that [A, [A,B]] = cil and [B, [A,B]\ = C2I for some 
ci,C2 € M, then the Baker-Campbell-Hausdorff formula^ reads 

(2.7) e^e^ = e^+^+II^-BJ+i^t^'I^.^ll-T^I^'I^.-B]]^ 
In particular, if [A, B] = cl then 

(2.8) e^+^ = e^e^e-^[^'^l 
Using (2.7) we have 

e-i?g(a+L)A ^ gL2(L+a)/6g(«+L)A+L(a+L)V-LC_ 

Using (2.7) again we deduce that 

g(a+L)Hg-L5g(a+L)A ^ ^L^a+L)/6 ^{a+L)H ^{a+L)A+L{a+L)V-Li 

By (2.8) we have e"«+(^'-"')^ = e"(^'-°')/2e°«e(^'-"')^, so the above identity gives 



Using this in (2.6) we deduce that 

Since e^^ is the shift operator {e'^'^ f){x) = f{x + r), we have e^^ Qm = Qrn-r/2^ ^"^^ we 
obtain 

The conclusion from the above is the following 
Lemma 2.1. 

We postpone the proof of this lemma until the end of this section. Putting this formula 
and (2.5) in (2.4) and using Lemma 3.1 of [CQR12] gives 
(2.9) 

P(sup(A(t) -t^) < m ) = lim P(^2(t) < for t £ [-L,a]) 

\t<a ) L^oo 

Having obtained an expression for the probability we are interested in, all that remains 
to show is that it coincides with our definition of (l-^)- We recall (as can be seen 

directly from its definition (1.9)) that the Airy kernel can be expressed as i^Ai = APqA* , 
where A is the Airy transform which acts on / G ^^(M) as 

/oo 
dzAiix-z)fiz). 
-00 

Since A* = A~^ we have by the cyclic property of determinants that the right hand side 
of (2.9) equals 

det(/ - PoA*Pj^+^2APo - Po^*e"^fe,+,2e-"fp^+„2^Po) • 



The Baker-Campbell-Hausdorff formula can be found in most introductory books on Lie groups and 
algebras. A general version can be found in [Dyn47]. However, in this very simple context, it is more readily 
computed by hand. 



and 
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Recalling that Qof{x) = f{—x) we have similarly that the last determinant equals 
det (^I - QopQ A* Pjri+a^APoQo - £'o^o^*e"^^+„2e""^P^+„2APo£'o^ 

Shifting the variables in the last determinant by —m we deduce that 

P( sup {A2{t) - t^) <m) = det(/ - PmEiPm - PmE2Pm) , 
\t>a J 

where 

dX M{x -m + 2m + 2a^ - A)e-2(^-^-°')" Ai(y - m + A) 

-00 

/>oo 

E2{x,y) = / dX Ai{x — m + X) Ai{y — m + X) . 

Shifting A by m + a^ in both integrals and changing A to —A shows that Ei{x, y) = K^{y, x) 
and E2 = K^, whence the equality in Theorem 1 follows since El = and E2 = K^- 

All we have left is to prove (2.5). We will denote by || • ||op, || • ||i and || • II2 respectively 
the operator, trace class and Hilbert-Schmidt norms of operators on L^(R) (see Section 3 
of [CQR12] for the definitions or [Sim05] for a complete treatment). Recall that 

(2.10) P^lli < PII2IISII2 and ||ylS||2 < p||2||S||op. 

Proof of (2.5). Let (p{x) = (1 + x"^)^^"^ and define the multiplication operator Mf{x) = 
ip{x)f{x). Then by (2.10) we have that 

Now ||e("+^)^/i:AiM-^||2 = ||M-ie("+^)^KAi||2 by the symmetry of e("+^)^ifAi, and then 
(3.3) in [CQR12] gives \\e^°'+^^" KaiM-'^^ < c{a + L)''^/'^ . Then to finish the proof it will 
be enough to estimate ||MP^^^2e~(°"''^^-'^-P^_,_(^2i^Ai||2 and \\MPjj^_^_j^2Ra^LPm+a''^Ai\\2- 

We start with the second norm. By (2.3), (2.10) and the fact that K^i is a projection 
we have 

\\MPrn+L'2Ra,LP—i+a'^KAi\\l < || MP^_,_2^2 i2Q^iP^_,_„2 || 1 1| i^Ai Hop 

Jm J -00 47r(a + L) 

where we have performed the change of variables x 1— >• x + L^. The inner Gaussian integral 
gives 

1 ^-4L3/3+2(«-L)x-4ma+4W+8aV3 N + ^rf/i + L)-'/'[Aa{a + L) + X - mlV 
4^7r(a + L) L ^2^ ^ ^ y ^ j;j 

where erf(z) = 27r-V2 f^dte'^'' < 1. Then 

/■oo 

||P^+^2i?„,,.P^+„2||2<C7e-4^^/3+4^"' /_ dx(^(x)2e2(°-^)-, 

•J m 

which clearly goes to as L — )• 00. 

On the other hand one can check that has integral kernel given by 

^-ia+L)H, . ^ 1 ^-Lx-a^+(a3+L3)/3-l^ 

^ V4vr(a + L) 

(this is done either by applying the Feynman-Kac and Cameron-Martin-Girsanov formulas 
as in [CQR12], or directly by integrating this kernel against the kernel of e("+^)^^ which 
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IS given, similarly to (2.1), by dX e-("+^)^ Ai(a; - A) Ai{y - A)). Note that this kernel 
is the same as the one given in (2.3) only without the reflection in the Gaussian term. It is 
easy to check then that the same calculation as the one in the above paragraph shows that 
||MP^+i2e-("+^)^P^+„2KAi||2 — ^ as L ^ 00. This finishes the proof of (2.5). □ 

Proof of Lemma 2.1. By (2.1) and (2.3), the kernel of e'^°^~^^^^ Kp^^Ra^L is given by 

/O roo 
d\ / dz e("+^)^ Ai(x - A) Ai(z - A) 
-00 J —00 

. ^ ^-{z+v-a'^-2rrf-L'^Y/A{a+L)+a^/i+L'i/i-Lz-ay 



V^47r(a + L) 

By completing the square in z in the exponential, the z integral can be seen as a heat 
kernel applied to an Airy function. Using the formula e*'^ Ai(x) = e^*' Ai(a; + t^) (see 
for instance Proposition 1.2 in [QR12]) we obtain after some manipulations 

i-G _ 

e("+^)^i^Ai^a,L(x, y)= d\ e2«'+2"'n^-2«?/ Ai(a; - A) Ai(2m2 ^2a^-y-\), 

J —00 

which corresponds to the claimed formula. □ 

We turn finally to the proof of Proposition 1.2. It relies on the following simple but very 
useful observation: 

Lemma 2.2. If A is a trace class operator on a Hilbert space %, 

det(/ + A) - 1 - tr(^) < ipilfeli^lii. 

Proof. Let A"'(^) = A® - ■ ■®A{n times, where ® denotes the tensor product of operators), 
which is an operator in the Hilbert space K^{T-L) known as the alternating product, see 
[Sim05] for more details. Then 

00 

+ A) = ^ti{K^{A)), 

k=0 

(this equality can be taken as the definition of the Fredholm determinant, see for instance 
(3.5) in [Sim05]). Of course A^{A) = I and A^(A) = A, so all we need to show is that 

00 

X^tr(A'=(^))<ip||?ell^lk 

fc=2 

But this follows directly from the inequality |tr(A^(74))| < ||A'^(^)||i < (/c!)^^ H^HJ' (see 
(3.1) and (3.4) in [Sim05]). □ 

Proof of Proposition 1.2. Using the inequality (Theorem 3.4 in [Sim05]) 

(2.11) |det(/ + A)- det(/ + B)\<\\A- ^H^e^+il^'li+H^ili 
for any two trace class operators A and we can write 

1 - Gl'^m) = 1 - det(/ - PmKlPm. - PmKlPm) 

< 1 - det(/ - PmKlPm) + ||P™/^2p^||^gl+2||P„i^ip„|U+||P„.i^2p„^||^_ 

We will show below that 

(2.12) \\PmKlPm\\i<ce-l'^''\ 

while a similar (and simpler) calculation gives the estimate ||-Pm-f^a-fm||i ^ ce~3'"^'' (this 
is in fact it is the same standard calculation based on (2.11) that gives the estimate 1 — 
-^gue(?7i) < ce~3™ ). Therefore all we need to show is that 

l-dei{I-P^KlP^)<ce- 



2 
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and thus, in view of Lemma 2.2, the proof it wih be enough to show that 

(2.13) tr(P^i^ip„) <cme-t'"'^'+-'" and \\PmKl^Pm\\i < ce-l^'^\ 
The trace can be computed directly: letting a = inax{0, a} we have 

triPmK^Pm) = dx I dAe°^Ai(x- A + a)Ai(x + A + a). 

Jm Jo 

We spht the integral according on the regions {A < x} and {A > x}. On the first region, 
changing variables A 1— >• xj, and since the Airy function is decreasing and positive on the 
positive axis and we may assume m > 0, the integral is bounded by 

dx (i7xe"^^Ai((l-7)x)Ai((l+7)x) < c / dx xe'"^T-3[(i-7f /^+(i+7)'/^]x^/ 

Jo Jm Jo 

where we used the first of the estimates 

(2.14) |Ai(x)| < ce"i'^''" for x > 0, |Ai(x)| < c for x < 
(see (10.4.59-60) in [AS64]). The term in brackets in the above estimate is larger than 2, 

and thus the whole integral is bounded by cme"'" 3" 
similarly, using (2.14), that the integral is bounded by 

/ dx (i7a;e"^^Ai((l -7)x + a) Ai((l + 7)x) < c / dx d-f xe'' 

J m J 1 J m J 1 

which again can be seen to be bounded by cme°'"~3™' . This gives the first bound in 
(2.13). 

For the second one write = {B'^e^^Po){Poe^^+°''>^B'^), where ^^(x. A) = Ai(x- A + a) 
and 5^ (A, y) = Ai{y + A + a). Then 

(2.15) WPmKlPmWi < WPmB^e-^PoMPoe^^+'^^^B^Pmh. 
The square of the first norm equals 

/"OO /"OO roo POO 

dx dXAi{x- X + afe^^^ < (ixe-2(^+"M dX Ai{- X)^ e"^^ < c e-"^"" , 

Jm Jo J m J — 00 

thanks to (2.14). Similarly, the square of the second norm equals 

dy dAAi(y + A + a)V(^+")^<c / dy e-i(^+^+")''^'+2(i+")^ < ce"!'"''^', 

J m J J m J 

again thanks to (2.14). Using these two bounds in (2.15) gives the second estimate in (2.13) 
and finishes the proof. □ 

Observe that for a < —1 the last bound in the proof above can be upgraded to 
ce^^ ' 3^ , givmg 



4 S/2 

and thus the whole integral is bounded by cme"'^~3^ . On the other region we have 

^ax7-|(l+7)3/2x3/2 



(2.16) WPmKP^t < ce"-"^-!'"''" > 0, 

which we used in justifying (1.10). A similar (although slightly more complicated) estimate 
gives ||Pm-^a-Pm||i ^ 0, which was used in the derivation of (1.11). 

a—)- 00 
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